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Abstract 

New inverse and half-inverse problems: sliding problems are intro- 
duced. In this way several physically important equations are recov- 
ered from the quantum defect. In particular, sliding problems are 
solved for radial Schrodinger equation, radial Dirac system and mul- 
tidimensional Schrodinger equation. Systems with Coulomb-type po- 
tentials are considered as well. 

In this paper we introduce new inverse and half-inverse problems: sliding 
problems. In Section [1] we show that this approach can be used for radial 
Schrodinger equation, radial Dirac system and multidimensional Schrodinger 
equation. In Section [2] we consider equations with Coulomb-type potentials. 

We note that inverse problems for the Dirac-type systems, which are 
usually called now Dirac systems, have been thoroughly investigated (see, 
e.g., [4^ITl [lTHl~4"ll20] and references therein). Radial Dirac system was intro- 
duced earlier than Dirac-type system and differs from Dirac-type system. 
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1 Inverse and half-inverse sliding problems 



1.1 Main definitions and results 

1. The radial Schrodinger equation has the form: 



dr 2 



+ 



+ 1) 



— q(r) ] y = 0, 0<r < oo (q = q), 



0,1,2, 



Some recent results and references on this well-known 



where * 

equation see, for instance, in [T||9|[TC)]. The radial Dirac system has the form: 

fi-(z + m-q(r))f 2 = 0, (1.2) 



d £ 
dr r 



d__l_ 

dr r 
0<r < oo 



f2 + (z- 

m > 0, 



m 



<l{r))h 
-q, £ = 



0. 



±1,±2,... 



;i.3) 

;i-4) 



where m takes positive (not necessarily integer) values and stands for mass. 
We easily rewrite (jl.2p and (11.31) in the matrix form 



d 
dr 



f = (i Z a 2 + V)f, 



1.5) 



/ = col [fx, fy], whereas the Pauli matrix a 2 and matrix function V have the 
form: 



0"2 



-i 

1 



V{r) 



—£/r m — q(r) 
m + q(r) Ijr 



;i.6) 



The peculiarity of the radial Dirac system consists in the fact that the entry 
ijr of the matrix V(r) is known and (£/r)^L 1 (0, oo). A well-posed inverse 
problem for the radial Dirac system was absent. This fact was discussed in 
our paper [T7]. In the present paper we introduce the notion of the quan- 
tum defect 8(r) at infinity (z = oo). The corresponding sliding half-inverse 
problem for Dirac system is well-posed. 

Without loss of generality we consider only the case i > since the 
equality (jl.5p yields the transformed equality, where Jf is substituted for 
/, — z is substituted for z, and — q and — £ are substituted for q and £, 

respectively (in the expression for V). Here J is given by J = J 
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Using quantum defect, we solve first a sliding inverse problem for the 
case of the radial Schrodinger equation. Consider equation (II .ip with the 
boundary conditions 

2/(0) = y(a) = 0, (1.7) 

and assume that 

/ \q(t)\dt < oo, < a < oo. (1.8) 
Jo 

The corresponding eigenvalues (i.e., the values of z such that the solutions 
y satisfying (\1.7\i exist) are denoted by z n (£) (n £ N). Then we have (as a 
corollary of Theorem 11.171) the asymptotic equality 

(L9) 

where 5(a), which is defined by the formula 

5(a) = [ q(t)dt, (1.10) 



is called the quantum defect of the equation (11.11) . Now, we formulate a 

sliding inverse problem: 

Problem 1.1 Recover the potential q from the given quantum defect 5. 

It is immediate from (11.101) that the solution of Problem 11.11 is given by the 
equality 

*M = QU) (.). (LID 

We deal with the radial Dirac system in a similar way. In order to for- 
mulate our result, we introduce a parameter e(z), which for the case of the 
Dirac system is given by e = \/m 2 — z 2 . 

Remark 1.2 More precisely, for all z G C excluding the cuts along the 



semi-axes (— oo, — m) and (m, oo), we choose a branch of e = ym" - '-' 



i 2 — z 2 

such that 31(e) > for z^O. We introduce e on the cuts as the product 
e = \Jm + z y/m — z, where 

, . i0, ifz>m, Jtt/2, ifz>m, 

arg(Vm + z) = \ arg(V«i -*)=<_ 

-7T/2, ij z < — m; 10, if z < —m. 



Hence, we have 

arg(e) = n/2, if z > m; arg(e:) = — vr/2, if z < —m. (1-12) 
The following theorem is proved in Subsection 11.31 
Theorem 1.3 Let the condition 

\q(t)\dt < oo, 0<r<oo (1.13) 

be fulfilled. Then there is a solution F(r,e,£) of system (jl.5p (where V is 
given by (11.61) ) . such that the relation 



F( r ,e,£)~--—=e-^ +s ^ 
2y/m + z 



e— >oo, z > m (1-14) 
holds and S(r) (i.e., the quantum defect of this system) has the form 

POO 

S(r) = / q(t)dt. (1.15) 



So, for Dirac system we formulate a sliding half-inverse problem: 

Problem 1.4 Recover the matrix function V in (11.51) from the given value 
I and quantum defect S. 

Since V has the form (ll.6p . the solution of Problem 11.41 is immediate from 
(11.151) or, more precisely, from the equality 

q{r) = ~ i'i 6 ) (r) - (L16) 

Remark 1.5 The notion of half-inverse problem and the first model of such 
problem were introduced in the paper fTh^ (see also JJE/). 

Remark 1.6 The values of 5(r) are spectral characteristics of the radial 
Dirac system. 

Remark 1.7 Sliding inverse Problem \l.l\ and half-inverse Problem \l.J\ are 

well-posed, that is, there is a one to one correspondence between the spectral 
data 5 and potentials q (see relations (jl.lOp and (I1.16p . respectively). Like 
all the classical inverse problems, these problems are not stable. 
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Remark 1.8 In the process of solving sliding Problems \1.1\ and 1.4 we use 
only the spectral data for positive energy. 



2. Next, we consider one- dimensional and multidimensional Schrodinger 
equations. Putting in ( 11. ip I = 0, we obtain the one-dimensional Schrodinger 
equation 

-y" + q{x)y = zy, 0<x<a, y{0) = y{a) = 0. (1.17) 

Theorem 1.9 ( [13j ) Let condition (11.81) be fulfilled. Then the eigenvalues 
Zn{a) = z n of system (11.171) satisfy the relation 

n 1 

yfz~=-n+- — 5(a) + o(n" 1 ), (1.18) 
a Znn 

where 5(a) is given in (ll.lOp . 

We note that formulas (II. 9p and (I1.18P coincide, when £ = 0. In Subsec- 
tion [L4] we prove the assertion: 

Theorem 1.10 Let condition (11.81) be fulfilled. Then the corresponding sta- 
tistical sum 

00 

Z{T) = Y,e~ Zn/T (1.19) 

n=l 

satisfies the relation 

Z(T) = (a/2) v / T/^- 1/2 + 5(a)/(2v / T^) +o(l/Vf), T^oo. (1.20) 

Recall that T stands in physical applications for temperature. 
The A;- dimensional Schrodinger equation has the form 

- Ay + q(x)y = zy, x = (x u x 2 , . . . , x k )eG C M k . (1.21) 

We deal here with a model case, where 

G = G{a u a 2 ,...) = {x : 0<Xi<ai (l<i<k)}, (1.22) 

q admits representation 

k 

q(x) = Y / q i {x i ) i (1.23) 
i=i 
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the functions % satisfy inequalities 

/ \qi(xi)\dxi < oo, < <2j < oo, l<«<fc, (1-24) 
Jo 

and the boundary condition (on the boundary T of G) is given by 

3/|r = 0, (1.25) 

Now, z„ j are the eigenvalues of the problem (I1.17p . where x = xi and q(x) = 
qi(xi), the statistical sum Zi is given by 

oo 

Zi(T) = ^V 2 "- l/T , (1.26) 

n=l 

are the points of the spectrum of problem (ll.2ip . (I1.25p . whereas the 
corresponding statistical sum Z(T) is defined again by the relation fll . 19j) . 
Representation (11.231) implies that 

k 

Z(T) = l[Z l (T). (1.27) 

i=i 

From Theorem 11.101 and equality (j 1.27ft we derive the following assertion. 
Theorem 1.11 Let the relations flL23l and (Ol hold. Then 

z{r > - {lY ( Vk - vf vi - + ^ + 2 ^ + ° ft)) • < L28 > 

where T— >oo and V_i = 0, Vo = 1, 

fc k 
V k = ]Ja i , V fc -i = Vfc^lK V k - 2 = V k Vffl^), (1-29) 

8=1 i=l l<i<j<fc 



5(G) = / g(x)da;. (1.30) 

./G 

The coefficients Vi have clear geometric interpretation, namely, 

Vi = Si/2 k -\ i = k, k - 1, k - 2, (1.31) 
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where Sk is the Lebesgue measure of the domain G, Sk-i is the Lebesgue 
measure of the boundary T, is the Lebesgue measure of the domain T 1 
formed by the set of pairwise intersections of the facets of the boundary T. 
Conjecture I. Formulas (ll.28l) - (ll.30l) hold for a much wider class of po- 
tentials than the potentials of the form (11.231) . 

Problem 1.12 Let formulas (ll.28l) - (ll.30l) hold. Recover the potential q from 
the given 5(G). 

According to (11.221) and (11.301) the solution of Problem 11.121 is given by 

Qk r 

q(a 1 ,a 2 , ■ ■ ■ ,a k ) = - — — / q(x)dx. (1.32) 



da\da 2 ■ ■ ■ dak 



G(ai,a2,...,ai c ) 



3. Finally, we consider one-dimensional and multidimensional equations 
for anharmonic oscillators. The one- dimensional equation for anharmonic 
oscillator (or one-dimensional anharmonic oscillator, which is shorter to say) 
has the form 

- y" + (x 2 /A + q(x)) y = zy, 0<x < oo (q = q), (1.33) 
and we introduce the boundary condition 

y(0)=0. (1.34) 

We assume that 

/■OO 

/ |g(a;)|(l + x 2+£ )dx < oo, e > 0. (1.35) 
Jo 

Theorem 1.13 ( [16J) Let condition (I1.35P be fulfilled. Then the eigenval- 
ues z n of problem f ll .33|) . ( 11.34j) satisfy the relation 

z n = In- 1/2 + —=(5 + o(l)), n^oo, (1.36) 
7tv 2n 



where 

poo 

q(x)dx. (1-37) 



Our next theorem easily follows from Theorem I1.13[ see Subsection 11.41 
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Theorem 1.14 Let an (inharmonic oscillator satisfy conditions of Theorem 
M.I'd Then the statistical sum Z(T), given by (I1.19p . satisfies the asymptotic 
relation 

Z(T) = T/2 -1/4-5 / (2VT^ + o (l/Vr\ . (1.38) 

The k- dimensional equation for anharmonic oscillator has the form 

f\x\ 2 \ k 

-Ay+ [LL. + g( V)J y = zy , xe G, \x\ 2 = £ X? . (1.39) 



We set 

G = {x : 0<Xi < oo (l<i<ib)}, (1.40) 

and introduce the condition (11. 25ft on the boundary T of the domain G. We 
assume again that the potential q(x) has the form ( 11.231) and that 

POO 

/ \qi(xi)\(l + Xi 2+£ )dxi < oo, e>0 (l<i<k). (1.41) 
Jo 

Like in the case of the multidimensional Schrodinger equation, z n ^ are the 
eigenvalues of the one-dimensional problem, this time of the problem (ll.33p . 
( I1.34p . where x = Xi and q(x) = q^Xi). By z n we denote the spectrum of the 
problem f l 1 . 3 9 [) with initial condition ( ll.25p . where G is given by (ll.40p . The 
statistical sum Zi(T) has the form (11.261) and the statistical sum Z(T) has 
the form (fl~T9]) . Like before, relation ( fl~23|) yields (IQTjl . From (Q7I> and 
(ll.38p follows the next assertion. 

Theorem 1.15 Let a k-dimensional anharmonic oscillator, where the poten- 
tial q admits representation (I1.23P and the functions qi satisfy the condition 
(ll.4ip . be given. Then we have 

-J (l -m/(2T)- 5(G) /(tV5v) +o(l/TVf)), (1.42) 

where T^-oo and 5(G) = j G q(x)dx. 

Conjecture II. Formula ( ll.42p holds for a much wider class of potentials 
than the potentials of the form (I1.23P . 

Remark 1.16 Theorem \1.15\ is an analogue of Theorem \l.ll[ Jointly with 
the results from ITSf . an easy modification of Theorem \1.15\ can be used for 
the formulation and solution of a sliding inverse problem for anharmonic 
oscillator. 
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1.2 Radial Schrodinger equation and quantum defect 

We add variables e = i\fz and I into the notation of the solutions of (11.11) 
and write y(r,e,£). Here, the parameter e, continuously depending on z, is 
defined for all z in the complex plane with the cut along the negative part of 
the imaginary axis. Let us fix arge = tt/2 for z > 0. We start our study of 
the equation (11.11) from the case q = 0. Two solutions ui(r, e, £) and u 2 (r, e, £) 
of (11.11) . where q = 0, are given by the relations 

Ul (r,e,£) = (2e)-^M . £+1/2 (2re), u 2 {r,e,l) = (2e) e W , i+1/2 (2re), 

(1-43) 

and M K ll (z) and W Kjfl (z) in the formula above are Whittaker functions. The 
connections between Whittaker functions and confluent hypergeometric func- 
tions are well-known (see, e.g., [5]): 

M 0>e+1/2 (x) = e-^V^a, c, x), W 0>e+1/2 (x) = e^V^a, c, x), 

(1.44) 

where a — £ + 1, c = 21 + 2. It follows from (11.441) and relations 

rfc i ) 

$(a,c, x)~l and ^(a,c, x)~ i=-. . x l ~ c for x-)-0, (1-45) 

r(a) 

that the solutions «i(r, e,£) and u 2 (r,e,£) satisfy the conditions 

Ul {r,ejW + \ u 2 {T,e,£)r I^±^- r-\ r^O. (1.46) 

Hence, in view of (I1.46P we have 

K(r,e,f)|<.M/ +1 for r<l/|e|. (1.47) 

We note that the letter Ai stands for different constants. 

Let us write the following integral representation (see [5j Section 6.11]): 

M 0/+1/2 (2re) = C(e)(2re) i+1 I(r,eJ), (1.48) 

where 

C(£) = 2-^)p|±||, I(r,e,£) = £^(l-tydt (1.49) 
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Integrating J(r, e, £) by parts £ + 2 times, we obtain the asymptotic formula: 
J(r, e, £)~{re)- {i+1) 2 l T{£ + 1) ((-l) m e - r£ 0(r, e, £) + e re 0(r, -e, £)), (1.50) 
where 

0(r,e,£) = l + ^ o +1) +o(l/e), £^oo. (1.51) 
ire 

Here, we used the Leibnitz differentiation formula in order to obtain an ex- 
pression for the {£ + l)'th derivative of the product (1 — t) e (l + t) e . Formulas 
(Q3D . fQ4l) and fO8D -f fT30l imply that 

Ml (r, e , £) = d (£, e) (-e irf / 2 e - r£ (r, e, £) + e^V^ (r, -e, £)) , (1.52) 

where d(£,e) = e i ^/ 2 ^g^(2 £ )-^+ 1 ). 

In order to study the asymptotic behavior of u 2 (r,e,£) we need the fol- 
lowing integral representation (see Section 6.11]): 

POO 

T(£ + l)W 0/+ i/ 2 (x) = e~ x/2 x e+1 / e~V(l + tfdt, x > 0. (1.53) 

Jo 

Formula ( 11. 53j) leads us to the equality 

POO 

T{£+ l)W 0/+1/2 (x) = e~ x/2 x- e / e- s s e (x + sYds, x > 0. (1.54) 
In view of the well-known integral representation 

POO 

T(z) = / e-V^ds, > 0, (1.55) 
</o 

equality ( ll.54p takes the form 

i 

T{£ + l)WW/ 2 (x) = e-^x^ ^ C p t Y{2£ -p + l)x p , (1.56) 

where C\ are the binomial coefficients. Using analyticity, we rewrite ( I1.56P 
for the complex plane 

i 

T(£ + l)W 0)t+1/2 (z) = e- z / 2 z~ e C ! T ( 2£ ~P+ l ) z ^ z ± °- (1-57) 

p=0 
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Now, the asymptotic relation 

u 2 (r,e,£) = (2e) e e- re (j)(r,e } e), e^oo (1.58) 

follows directly from fl 1 . 43 1) and (ll.57p . 

Next we deal with the case q^O. We assume that q satisfies the condition 

/ \q(t)\dt < oo, 0<r<cx). (1.59) 
Jo 

Then there is a solution u of (II. ip such that the equality 

u(r, e, £) = «i(r, e,£) + / G(r, t, e, £)q{t)u{t, e, £)dt (1.60) 

Jo 

holds for G of the form 

G(r, t, e, i) = {Mr, e, £)u 2 (t, e, I) - u 2 (r, e, £) Ul (t, e, £)) . (1.61) 

Relations ( I1.46P and fll.6ip imply that 

{t/r) i+1 \G{r,t,e,£)\<Mt, t<r<l/\e\. (1.62) 

From fll.46p . (I1.60p and (ll.62p we derive that the asymptotics of u is similar 
to the asymptotics of U\, namely 

u(r,e,£)~r e+l , r<l/|e|. (1.63) 

Moreover, according to (ll.5ip . (I1.52p . (11.581) and (11.611) we have 

G(r,t,e,£)~— (e (r -* )e -e- (r - t)e ), s^oo. (1.64) 

Using (I1.46P and (ll.60p - (ll.63p we obtain the relation 

•i/kl 

\G(r,t,e,£)q(t)u(t,e,£)\dt = o(\e\- (e+2) ) , e^oo. (1.65) 

It follows from (ll.64p that for some Ai we have 

\G{r,t,e,£)\<Me rlal /\e\, l/|e|<f<r, (1.66) 
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o 



where a := 8(e). In view of (OS) . (fl~5T) . Q2} ; (OS]) and (Q3j) the 
relation 

1/1=1 

\G(r,t,e,£)q(t)u(t,e,£)\dt = e r|<T| o(l/|£|), r > l/\e\, e-^oo (1.67) 

holds. Formulas (11.60) . ( jl.66p and (11.67) imply the following equality 

w(r,e,£) = Mi(r,e,£) + (e rH /N) , r > l/\e\. (1.68) 

We denote by z n (a, £) the eigenvalues of problem ( 11. ip . (11.7) . It is well-known 
that there is only one (up to a constant factor) solution of ( II. ip which turns 
to zero at r = 0. Since u is such a solution, z n (a,£) is an eigenvalue if and 

only if 

u(a, iy/z n (a,£), £) = 0. (1.69) 

By well-known methods (see, e.g., [121 Ch.l, Section 2]), using relations 
( E2D , (USED and (OSl . (Q5)l . we obtain that 

V^M) = ^ (™ + + o(l). (1.70) 

Recall that a = 9ft(e). Putting 

»(e) = 0, (2e) m M (r, e, £) = u(r, e, £), (2eY +1 Ul (r, e, £) = u x (r, e, £), 

and taking into account estimates ( 11.64) and ( 11.681) . we deduce from ( 11.60) 
that 

u(r, e, £) = ui(r, ej) + — [ (e (r "* )e - e~ ir - t)£ ) ut(t, e, £)q(t)dt + o(- 

2 £ Jo \ £ 

(1.71) 

Thus, we obtain the following assertion. 

Theorem 1.17 Let condition (11.591) be fulfilled and equality 3ft(e) = hold. 
Then the solution u(r,e,£) of the equation ( 11.1) satisfies ( 11.631) and admits 
representation 

U(r,e,£) = 0,(1,6) (_ e ^/2 e — +A(r)/(2 £ ) + e - W /2 erc -A(r)/( 2£ ) + (l/| £ |)) ? 

, (1-72) 

where A(r) = -£(£ + 1) — 8(r) and the quantum defect 8(r) has the form 

rr 

8(r) = / q(t)dt. (1.73) 
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Recall that ( II. 9p follows from Theorem II. 171 We denote by z n (£) the positive 
roots of the Whittaker function Mo^+i/2(i^), which are not equal to zero. 
Formula (II. 9p yields the next well-known result (see [HI Ch. XI]): 

^ H .(„ + 0_J±i2_ + (1.74, 

1.3 Dirac equation and quantum defect 

1. We assume that q satisfies (jl.4j) and (jl,13p . that is, 

/"OO 

q(r) = q(r), / \q(t)\dt < oo, r > 0. (1.75) 

Given functions /i and f 2 , it is easy to find functions Q\ and Q2 such that 

fi(r,e,£) = V^T^e- r£ (2re) e - 1 (Q 1 (r,e } £) + Q 2 (r,e,£))r, (1.76) 
f 2 (r,e,£) = -V^^e- r£ (2rey- 1 (Q 1 (r,e,£)-Q 2 (r,e,£))r, (1.77) 

where e = y/ m 2 — z 2 . Recall that e is defined more precisely in Remark 11.21 
In order to prove Theorem 11.31 following [2J we consider solutions of (II. 2\\ - 
(OD in the form (TL76|) and (TLTTI) . Substituting (TL76|) and (TIT?]) into (TOj) 
and (11.31) we obtain: 



/ TOO V 

r(Qi + Q 2 )' + 2£(Q, + Q 2 ) - 2rsQ 2 + J — (Q 1 - Q 2 )rq(r) = 0, (1.78) 

m + z 



I rji _)_ y 

r(Qi - Q 2 )' + 2reQ 2 - J -—- {Q 1 + Q 2 )rq(r) = 0, (1.79) 

V m — z 

where Q' = 4-Q. It is easy to see that 



m — z m + z 2 m — z m + z 1 

J = — z, J + J = -m. 1.80 

m + z \ m — z e \ m + z y m ~~ z £ 

Hence, if q(r) = 0, the equations ( I1.78P and (I1.79P can be rewritten in the 
form 

rQ[ + £Q X + £Q 2 = 0, rQ' 2 + £Q 2 - 2reQ 2 + £Q X = 0. (1.81) 
We substitute into ( ll.8ip functions 

Q l (r,e,£) = Q h0 (p/(2e),ej) (i = 1,2), p = 2re, (1.82) 
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and so switch to the new variable p. It follows from (jl.8ip that 

Pj^Qifi + (2* + 1 - P)j-Qifi ~ *Qi,o = 0, (1.83) 

pjfQw + ( 2i + 1 - p)j-Qw ~ ( £ + VQw = °- ( L84 ) 

The solutions of (I1.83P and ( ll.84p . which are regular at the point p = 0, are 
the confluent hypergeometric functions ai$(£, 2£ + l,p) and a 2 §(£ + 1, 2£ + 
l,p), respectively (see [2] Section 36] and [3, Section 14]). Substitute r = 
into (jl.Bip . f ll.82p to see that a\ = —a 2 . That is, up to the factor ai we have 

Q lfi (p) = m 2£ + 1, p) = e^ 2 p- e - 1 / 2 M 1/2ie (p), (1.85) 
Q2M = ~m + 1, 2£ + 1, p) = -e^ 2 p^ 1/2 M_ 1/2/ (p). (1-86) 

Now, we need an integral representation from [5J Section 6.11]: 

M ±1/2 ,(p) = C(£)(pY +1 / 2 I ± (p,£), C(e) = 2- e T T W+V (1.87) 
where 

I + {p,£) = ^ e tp/2 (l-t 2 y-\l-t)dt, I_{p,£) = ^ e tp/2 (l-t 2 ) e -\l+t)dt. 

Integrating the expressions above by parts (£ times), we obtain: 

I + (p,£M-lY(p/2)- £ T(£)e- p / 2 , p^oo (arg(p) = vr/2), (1.88) 
/_(p,£)~-(p/2)^r(£)e"/ 2 , p^oo (arg(p) = vr/2), (1.89) 

Using asymptotic formulas ( I1.88P and (I1.89P we have 

Q lj0 (p) ~ (-l)V'r(2£+l)/r(^+l), p^oo (arg(p) = tt/2X,1.90) 
Q2M p- e e p T{2£+l)/T{£+l), p^oo (arg(p) = vr/2). (1.91) 

Formulas ffL90|) and f ll~9Tj) show that the solutions of ffL83|) and ffl~Mj) of 
the form (I1.85P and (I1.86P are, indeed, regular at p = and, moreover, we 
have 

<2i,o~l, Q 2 ,o~l for p^O. (1.92) 
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The functions 

Q 1)0 (p) = Cl *(e, 2i + 1, p) = dePPp-'-^Wx/zM (1-93) 
Q2M = ^{£ + 1,2£+ 1, p) = czer^p-'-^W-wM, (1.94) 

where ^(a, c, x) is the confluent hypergeometric function of the second kind, 
also satisfy (11.831) and (11.841) but are non-regular at p = (see [5], Ch.6]). 
According to [5l Section 6.11], we have the following integral representation 
for W ±1 / 2 ,t- 

W ±l/2/ {p) = C ± {i)e~p' 2 {pY +1 / 2 J±(p, I), C ± (£) = 1/T(£ + 1/2 T 1/2), 

(1.95) 

where 

poo poo 

J + (p,£) = e- tp t e -\l + t) e dt, J_(p,£)= e- tp t\l + ty- l dt. (1.96) 
Jo Jo 

Using integral representation ( ll.55p . we rewrite (I1.96P as 

i i-\ 
J + (p, £) = p- 21 QT(2£ - i)p\ J_(p, £) = p~ u £ Ci_ x T{2l - i)p\ 

i=0 i=0 

(1.97) 

where C\ are the binomial coefficients again. The right-hand sides in (I1.97P 
are analytic functions (for p^O). Thus (11.971) is valid for all p^O in the 
complex plane. 

In view of (Q5) - (Q51l and fTOTj) . if we have 

QiM-c lp - 2e T(2£)/T(£), Q 2>0 (p)~ C2 p- 2e r(2£)/r(£ + 1). (1.98) 

After the substitution Qi,o = Qi,o an d Q2,o = Qi,o into (11.821) . relations 
(ll.8ip and (11 .98[) imply that we can assume 

d = 1, c 2 = I (1.99) 

Using again formulas (ll.93l) - (ll.95l) and (ll.97p . we derive 

g li0 (p)~p- £ , Q 2 ,o(p)^P~ (m) , P^oo (arg(p) = 7r/2). (1.100) 

The constructed regular and non-regular (at r = 0) solutions of (ll.2l) - (ll.3l) 
(or, equivalently (11.51) ). where q = 0, we denote by 

F (r) = F {r,£,£) = col [F 1>0 (r, e, £) F 2)0 {r,£,£)\ (1.101) 
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and 



(1.102) 

(1.103) 
(1.104) 
(1.105) 
(1.106) 

mv+i) (L107) 

xcol [((-l)V r£ -e r£ )V^T^ - ((-l) e e- r£ + e r£ )Vm- z\ , 

where z > m (i.e., arg(e(z)) = vr/2) and either e — > oo or r — > oo. In a 
similar way, from ( 11. 1001) we derive that 



F (r) = F (r,M) = col [F 1)0 (r, e, F 2>0 (r, e, £)}, 

respectively. According to (I1.76P and (ll.82p we have 

*i, (r) = v / ^T^e- r£ (2re) f - 1 (Q li o(r) + g 2 , (r))r, 
F 2 , (r) = -V^^e- r£ (2reY- 1 (Q lfi (r) -Q 2fi (r))r; 
F 1)0 (r) = v / ^T^e- r£ (2r£) f - 1 (Q 1 , (r) + Q 2 , (r))r, 
^a,o(r) = -v / ^^e- r£ (2re)^ 1 (g li0 (r) - Q 2 ,„(r))r. 

Hence, from (jl.90p and (ll.9ip we obtain 



F„(r,£,f)~-e- re col 



m, 



(1.108) 



where either e — > oo or r — >■ oo 

2. Now, we consider the case when q(r) ^ 0, and introduce 2x2 matrix 
functions 



V(r) 



q(r) 
—q(r) 



, U (r,e,. 



F 1>0 (r,e,£) F li0 (r,e, 
F 2i o(r,e,£) F 2 , (r,e, 



It is easy to see that the solution F(r,e,£) of the equation 

/■oo 

F(r) = F (r) + / U^U^V^F^dt 



1.109) 



(1.110) 



satisfies the differential system (ll.2p . (jl .3j) . 
The following equality is valid: 



dettfofoM) = r(2£)/( £ r^)). 
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(1.111) 



Indeed, it follows from (ll.2p and (jl.3p that ( det U (r, e, £))' = 0. Hence, we 
obtain the relation det Uo(r, e, £) = const. Since the first and second columns 
of Uq are F and F , respectively, formulas fl 1 . 107|) ( ll . 109j) imply that 



det U (r , e , t) ~r (2£) / (eT (£) ) , r— >oo (z > m), 

which yields (jl.llip . 

According to (O07ll - (ll~T09ll and (OTTl) . the equality 

f/o(r)t/ (t) _1 = e< r -*> e ei + e (t - r)e 6 2 + o(l), £^00 (2 > 



m 



1 i 
-i 1 



1 -i 
i 1 



;i.ii2) 
;i.ii3) 

(1.114) 



holds. We multiply both sides of ( 11. 11011 by -2y/m + ze re . Using ( 11. 113ft 
and passing to the limit z— >■ + oo, we obtain 



i 








1 





where 



The equality 



—2 lim ( \/m + ze r£ F(r, e, 



-if r °°<l(t)dt 



;i.ii5) 
;i.ii6) 

1.117) 



follows directly from (ODBjl . (005)1 . (OI3)l and ( TTTTBll . From flTTTKH 
and ( II. 117ft we see that the constructed solution F(r, e, £) has the properties 
stated in Theorem 11.31 that is, Theorem 11.31 is proved. 



1.4 Proofs of Theorems 11.101 and 11.141 

In order to show that the results on multidimensional Schrodinger equation 
from Subsection 11.11 hold, we should prove Theorems 11.101 and 11.141 First, 
we prove the assertion below (see [T9||2T]). 

Lemma 1.18 The following asymptotic relation 

oo 

^ e -^ = -v^F-- + o(l), z^oo (1.118) 

71=1 

is valid. 



17 



Proof. In the proof of (11.1181) we use the Poisson formula (see, e.g., [B]): 

00 1 POO 00 POO 

Ve"" 2/z = — +/ e- x2/z dx + 2j2 e- x2/z cos(2irnx)dx. (1.119) 
2 J Q J 



n=l 

Since 



POO 1 

/ e- x2 / z cos(27mx)dx= -^e~ zn27T \ (1.120) 
Jo 2 

formula fl 1 . 11 9 [) implies that 



J>-" 2 /* = -- + -v^F+ v^E e " 2nV - ( L121 ) 

n=l n=l 

Relation (jl.lisp follows directly from fll.mp . □ 
Proof of Theorem \1.1(A In view of f 1 1 . 1 8 j) we have 

Zn = (nir/a) 2 + (l/a)5(a) +o(l), n-kxx (1.122) 

Substituting z = T(a/7r) 2 into (1 1 . 1 1 8 [) and taking into account (11.1221) we de- 
rive (II . 20 j) . where Z is given by f II . 1 9 j) . □ 
In the proof of Theorem I1.14[ the next well-known relations 

i i 



S e " 2 " /T =2imh(l7T) = 2 (r + 0(1/T »- T ^ < L123) 

n=l v ' 7 

OO /»CXJ 

Ve- 2n/T /72^= / (e- 2a; / T /V2^)dx + 0(l), T^oo (1.124) 
n=l ■'o 

are used instead of the Poisson formula, which was applied in the previous 
proof. Since 

POO 

/ ((T u /\fu) du = v 7 ^, (1.125) 
Jo 

the assertion of the Theorem 11.141 follows from (I1.123P and (I1.124p . 
1.5 Dirac system on a finite interval 

The considerations from the previous subsections can be used also to deal 
with Dirac system on a finite interval. In this final subsection of the section 
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we consider the Dirac system (11.21) . (11.31) on the segment [0, a] and assume 
that the inequality 

/ \q(t)\dt <oo (1.126) 
Jo 

holds. The solution F(r, e, £) of the integral equation 

F(r,e,£)=F (r,e,e)- [ U Q {r)U Q {ty l V{t)F{t,eA)dt, (1.127) 

Jo 

where Uq and V are given by (11.1091) . satisfies the system (11.21) . (II. 3p . Here 
F is a regular solution of the system (II. 2p . (II. 3p with q = and is given by 
(11.1011) . (11.1031) and (II . 104j) . The parameter e = e(z) above is the same as in 
Remark 11.21 and in Subsection 11.31 Recall that -8(0, I] stands for the class of 
bounded functions on (0, I], that is, such functions, the values of which are 
uniformly bounded in the norm. Below we consider functions bounded with 
respect to the variable r. From formulas (11.851) . (I1.86j) . (1 1 . 1 3 j) and ( II . 104j) 
we derive 

F (r,e,£)/ ( v / ^T^(2re)^ 1 r) e B(0, l/\e\] (z > m). (1.128) 

In fact the left-hand side of (1 1 . 128j) is bounded when z > m and p = 2re is 
bounded. In a similar way, using (11 .93j> . ( j 1.94ft . (II. 105ft and (11.106ft . we have 

F (r,£,£)/ (V^Tz(2re)- e - 1 r) e 5(0, l/\e\) (z > m). (1.129) 

According to (ODD . (OTj) and (fi~98l) the relation 

U (r) E B(0, l/\e\] (z > m) (1.130) 

is valid. Relations (fl~T09j) . (Tl~T28|) . (029]) and (Tl~T30|) imply that 

^T^{2rey- 1 rU a {r)- 1 e B(0, (z > m) (1.131) 

and 

sup||(t/r) 2 ^ (^)^oW _1 || < oo, 0<t<r<l/|e| (z>m). (1.132) 
It follows from (lITSTj) . (028]) and (0321) that 

F(r,e,£)/ ( v / ^T^(2re)^ 1 r) E [B(0, l/\e] {z > m). (1.133) 
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In view of (O30) . (Q3I]) and (035|1 the relation 

■i/kl 



Uo(r)U (t)- 1 V(t)F(t,£,e)dt 



o(l), £ -)■ 00 (2 > m) (1.134) 



is valid. Substituting fTOOj) and f TOTj) into (jl.l03j> and fll.l04jt . we represent 
i*o(r, £,£) (for the case that e = e(z), z > m) in the form 



oir,£, 



e -V £7r/ V +e r£ e- i£7r/2 ^ 2 + o(l), £^00, (1.135) 



£1 :r 



r(2£ + i) . 



2T(£+ 1) 



92 := rg+^e^ 



2r(£ + 1) 



(1.136) 



Now, we can prove our next theorem. 



Theorem 1.19 Let condition (I1.126P hold. Then (for the case that e = e(z), 
z > m) we have 



e-^ 2 F(r,£, 



e" r£ e 



xK) - 2r(*+i) 



ife/2 /ii(r) + e r£ e-^ /2 /i 2 (r) + o(l), £ ->• 00, (1.137) 

i 



2r(£ + 1; 



where the quantum defect S(r)has the form 



6(r) 



q(t)dt. 



1.138) 



;i.l39) 



Proof. In order to prove the theorem, we represent F(r,e,£) in the form 

e -^/ 2 F(r,e,£) = e - re e i£7r/2 /n(r) + e r£ e^' 2 h 2 (r) + f(r,e,£) (1.140) 
and estimate /. Definitions (11.1361) and (11.1381) imply that 

h 1 (r) = g 1 - f ' TtVtytuWdt, h 2 (r)=g 2 - f T 2 V{t)h 2 {t)dt. (1.141) 
Jo Jo 

Let us multiply both sides of (11.1401) by ^l 2 and substitute the result into 
(11.1271) . Then, formula (11.1411) implies that the function 

G(r, £, £) = f(r, e, £) + f U^U^V^f^ e, £)dt (1.142) 

Jo 
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satisfies the relation 

\\G(r,e,£)\\->0, z^ + oo. (1.143) 
Using ffl~T26|) . (1TT42]) and (OI3]) . we obtain 

||/(r,e,£)H0, z^ + oo. (1.144) 

□ 

Recall that F = col [Fi F 2 ~\ is the regular solution of the radial Dirac 
system. Consider the case of the second boundary condition 

F(a, e n , £) sinip + F 2 (a,e n J) cosip = 0, -vr/2<^<7r/2. (1.145) 

Here e n = e(z n ) = a/ m 2 — z\ and (differently from other considerations of 
this section, where n always belongs N) n G Z — {0}. Without loss of 
generality we assume that z n > m for n > and z n < —m for n < 0. 

Corollary 1.20 Lei conditions (11.1261) and ( 11. 145ft be fulfilled. Then we 
have 

7T / A V - + 5(a) , , , 

z n (a) = -[n+-)+- — + o(l). (1.146) 



Proof. It follows from Theorem 11.191 that 

i^r,^) = ~e ifa/2 ^^ sin (rs n - fa/2 - o(r)) + (1), (1.147) 

F 2 (r,e n ,£) = -e^ /2 r ( (^+i) ) cos ( rs " " £7r / 2 " + ( L148 ) 

where e n = iW z\ — m 2 = is n (n > 0, z n > m, s n > 0). In view of f l 1 . 1 45 j) . 
(fj~T4"TD and (frA . the relation 



cos 



(as n - £tt/2 - 6(a) - iji) = o(l) (1.149) 



is valid. Hence the equality ( I1.146P is proved for n > 0. The case n < can 
be dealt with in the same way. □ 

Remark 1.21 Formula ( 11.1461) is essential for solving the corresponding in- 
verse sliding problem. 

Remark 1.22 When £ = 0, formula (11.1461) is well-known (see JT^ Ch. VII]). 
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2 Schrodinger and Dirac equations 
with Coulomb-type potentials 

The radial Schrodinger equation with the Coulomb-type potential has the 
form 



(Py ( 2a £(£+1 
dr 2 V t r 



2 + U + -^-q{r))y = 0, 0<r<oo, a = a ^ 0, (2.1) 



where £ = 0, 1, 2, . . . Its potential differs from the potential in (II. ip by the 
additional term — . The radial Dirac system (relativistic case) with the 
Coulomb-type potential has the form: 



{t + ^) h ~ (^ + ™ + ^-<?(0)/ 2 = o, 



(2.2) 



'' ^ f2+(z-m + --q(r))f 1 = 0, (2.3) 



dr r J V r 

0<r < oo, m > 0, q = q, a = a^0, £ = ±1,±2,..., £ 2 >a 2 . 

(2.4) 

Equations ( 12. 2p and < \2.3h differ from the equations ( 11. 2ft and ( 11. 3ft by the 
additional term - . Like in Section HJ without loss of generality we consider 
only the case £ > 0. 

The scheme presented in Section [1] admits modification for the case of 
Coulomb-type potentials. Equation (12. ip and system ( 12.11) . ( 12. 2 p are essential 
in the study of the spectrum of atoms and molecules [21131122] . 

We assume that 

\q(t)\dt < oo, r>0. (2.5) 

In the present section we describe the asymptotic behavior of the solutions 
of equation (12. ip and of system (I2.2p , (12. 3p for the energy tending to infinity 
(i.e., z— >oo). Using this asymptotics, we introduce the notion of the quantum 
defect and solve the corresponding sliding inverse problem for the relativistic 
case. 
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2.1 Asymptotics of the solutions: Schrodinger equa- 
tion 

If q(r) = 0, then ( 12. ip takes the form: 

d 2 y ( 2a £(£ + 1) 



dr 2 \ r r 2 j ^ ^ 

We denote by Ui and u 2 the solutions of (12.61) such that 

«i(r, e, i) = (2s)- (e+1) M a/e/+1/2 {2re), u 2 {r, e, £) = (2eYW a/£/+1/2 {2re), 

(2.7) 

where M Ktfl (z) and W K)ll (z) are Whittaker functions. Like in Subsection II. 2\ 
the parameter e = iy/z is defined for all z in the complex plane with the cut 
along the negative part of the imaginary axis. We put arge = n/2 for z > 0. 
Recall the connections between the Whittaker and confluent hypergeometric 
functions (see [5]): 

M a/£/+1/2 (x) = e- x/2 x c/2 $(a, c, x), W a/e/+1/2 (x) = e~ x/2 x c/2 ^(a, c, x), 

(2.8) 

where a = £ + 1 — a/e, c = 2£ + 2. Instead of asymptotics in (11.45j) . we 
need now the asymptotics of $ and \I> for energy tending to infinity. More 
precisely, we need the following well-known formulas (see [5j Section 6.13]): 

*(a, c, 2re)J-^±^-{-2re)-^ + J^^^e)-^e 2 -, (2.9) 

'^(a,c,2re)^(2re)^^ 1 , where e— >oo and aig(e(z)) = n/2 (i.e., z > 0). 

(2.10) 

Using (12. 9p and ( I2.10p . we obtain the relations 

Ulir ' £) ~ T T(e + l) ((- 1 ) me " r£ + e re ) (2^)- (m) , e->oo, z > 0, (2.11) 
u 2 (r,e)^(2£Ye- r£ , e^oo, z > 0. (2.12) 

Next, we consider the general- type Schrodinger equation ( 12. ip (where q 
is not necessarily trivial). The solution u(r,e,£) of the integral equation 



u(r, s, £) = u 2 (r, e,£)- / k(r, t, e, £)q(t)u(t, e, £)dt, (2.13) 

r 
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where the kernel k(r,t,e,£) is defined by 

k(r, t, e, I) = r ^| 2 £ + 2) /£ ^ £ ' ^ U2<(t ' £ ' ^ _ u ^ r ' £ ' ^ Ul (*> £ ' ^)) ' 

(2.14) 

satisfies fl2ZE|). It follows from ( 12TTT) . (12TT2D and (IZHjl that 

fc(r,t,e,£)~e (r - t)£ -e- (r -* )e , e^oo, 2 > 0. (2.15) 

Using standard methods, we deduce from (I2.12I) - (I2.15I) the following state- 
ment. 

Theorem 2.1 Let (12. 5p /ioW. Then, for u(r,e,£) satisfying (I2.13P we /iai>e 
M(r,e,£)~(2e)-V re , e^oo, z > 0. (2.16) 

2.2 Asymptotics of the solutions: Dirac system 

If q{r) = 0, system §F2$, (£3} takes the form 

f d £\ „ . a. „ ( d £\ „ , a. „ ,„ 9 9 . 

(* + r) = "-(*-;) = > «')■ 

(2.17) 

We consider solutions of (I2.17P in the form similar to (I1.76P and f ll.77p . that 
is, 

A(r, e, £) = ^n~T~z e^r^- 1 (Q 1 (r, e, £) + Q 2 (r, e, £))r, (2.18) 
f 2 (r, e, £) = -V^~ze- re {2reY- 1 (Q 1 (r, e, £) - Q 2 (r, e, £))r, (2.19) 

where u = \/ £ 2 — a 2 > 0, e = \fm 2 — z 2 and the choice of arguments in 
y/m ± z and \/m 2 — z 2 is prescribed in Remark 11.21 For regular and non- 
regular solutions of (I2.17P we use here the same notations as for regular and 
non- regular solutions of (11.21) . fll.3|) . where q = 0, in Section [TJ For the case 
of regular (at r = 0) solutions of (I2.17p . the functions Q\ and Q 2 can be 
expressed via the confluent hypergeometric functions $(a,c, x) (see [2J[3]). 
Thus, a regular solution F of (I2.17P is given by 

F = col [F lfl F 2)0 ] , (2.20) 
F 1)0 (r, e, I) = V^+~z e- r£ (2rer- 1 (Q 1)0 (r, e, I) + Q 2>0 (r, e, £))r, (2.21) 
F 2 ,o(r,e,£) = -V^^e^re)"- 1 ^^, M) - Q 2>0 (r,e,£))r, (2.22) 
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where 

Qi,o = ai&(u — az/e, 2u + 1, 2re), Q 2: o = « 2 $(w + 1 — az/e, 2u + 1, 2re), 

(2.23) 

and 



?1±^ = - ° /l^. (2.24) 

Using asymptotic formulas (see [5]) for the confluent hypergeometric func- 
tions $(a, c, x), we obtain 

Q 1)0 (r, £,£)-«! -I^±lL(2 re )--- ia , e^oc (z > m), (2.25) 
r(a; + 1 — la) 

Q 2 , (r, £ ,£)^a 2 r . r(2 ^ i + | 1) . e 2r£ (2r g )- +iQ , e^oc (z>m), (2.26) 
f (a; + 1 + ia) 

where 

g 1 + g 2 = _-2!L_. (2.27) 
«i — a 2 w + £ 

A non-regular at r = solution F = col F 2t o] of ( I2.17p has the form 
(see 0E]): 



F lj0 (r,M) = V^T^e- re (2rer- 1 (Q lfi (r,e,e) + Q 2j0 (r, M))r, (2.28) 
^2,o(r,£^) = -V^^e- r£ (2rer- 1 (Q lfi (r,eJ) - Q 2fi (r,e,£))r, (2.29) 

where 

Q 1)0 (r, e, £) = 5 x *(w - az/e, 2a; + 1, 2re), (2.30) 
g 2 , (r, e, £) = a 2 ^(co + 1 - az/e, 2tu + 1, 2re). (2.31) 

Here ty{a,c,x) is the confluent hypergeometric function of the second kind. 
In view of fl2U0|) and (Oljl we have (see Ch.6]): 

Qi,o~Si(2re)~ 2a; , Q 2 ~5 2 (2re) _2w , r -> oo; 5?i = — — —a 2 . (2.32) 

we + az 

Using again asymptotic formulas for c, x) (see |5]), we obtain 

Qi,o~«i(2re)" a; ~ ia , e— t-oo (z > m), (2.33) 
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Q 2 ^a 2 {2re)-^ a -\ e^oo {z > m), (2.34) 

where 

t u . 

«i = «2- (2.35) 

u — la 

From ( l2T2Tj) . fl2T22|) . f l2T25|) and fl2T26|) we obtain 



F li0 (r, e, £)~VrrT+z T{2u + 1) 



x j e — a 1 f 2r£) - . + e r£ a 2 j^^- T ) . (2.36) 

1 2eT(u + l-ia) 2eT(cj + 1 + id) 1 



F 2 , (r, e, e)~ - T(2w + 1) 

/ (2re)- ia re ^ (2re) ia \ /n _ 

x e re ai— -j - e rE a 2 —^r z r , (2.37) 

V 2eT(u + l- ia) 2eT{u + 1 + ia) J ' v ' 

where e^oo, z > m. Formulas f l2T28|) . (I2T29D and (T2T331 - (12351) imply that 

F (r, e, £)~—aie~ r£ (2re)~ ia col [y/m + z — y/m — z], e— t-oo (z>m). 

(2.38) 

Next, we consider the case q(r) ^ 0. The solution F(r, e, £) of the integral 
equation 

POO 

F(r 1 e,£)=F (r 1 e,£)+ / U Q (r)U*{t)- x V{t)F{t,e,C)&, (2.39) 



where V and Uq have the form ( 11 . 109H . satisfies system (12. 21) . ( 12. 3ft . We note 
that -F^o and F^o (i = 1,2) in Section Q] are different from the entries F^q 
and Fj of C/q, which are introduced in this section. For the present case, 
formulas (12361) - ( 12381) imply that 

det£/o(r,g,l)~-a 2 Q!i n X t • v e_K3 ° ( z > m )- ( 2 -40) 

2el (a; + 1 + ia) 

According to (12T361) - ( 12381) and lj£30jl the equality IjTXTSjl from Section Q] is 
valid again. Hence, multiplying both sides of (I2.39j) by — 2-y/m + i e r£ (2re) m 
and passing to the limit z— > + oo, we obtain 



Foo(r,. 



i 




1 


+ / 



+ / Q.V^F^^dt, (2.41) 
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where 

F oc (r,£) = -2 lim (V^h~T^e re (2re) ia F(r,eJ)). (2.42) 

z— >+co \ / 

The equality 



F 0O (r,£) = e- i f?' q ® dt 



(2.43) 



follows directly from ( I2.4ip -( )2"T4"3"]) . Thus, we proved that F(r,e,£) con- 
structed above satisfies the requirements of the following statement. 

Theorem 2.2 Let condition ( 12. 5 p be fulfilled. Then there exists a solution 
F(r,e,£) of system (12. 2p . (12. 3p . which satisfies the relation 



F(r, e,£) jl e - i(re+5(r)) (2re) 

2\/m + z 



, e->-oo (2; > m), (2.44) 
where the quantum defect 5(r) is given by the formula 

POO 

S(r) = / q(t)dt. (2.45) 



Finally, we formulate the sliding half-inverse problem. 

Problem 2.3 Recover the potential q — a/r of the Dirac system (12.21) . ( 12.31) 

from the given quantum defect 5 and constant a. 

According to (12.451) . the solution of Problem 12.31 has the form: 

q {r) = ~^8{r). (2.46) 

til 
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